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Mediterranean climate regimes are particularly susceptible to rapid shifts between drought and flood—of which, California's 
rapid transition from record multi-year dryness between 2012 and 2016 to extreme wetness during the 2016-2017 winter pro¬ 
vides a dramatic example. Projected future changes in such dry-to-wet events, however, remain inadequately quantified, which 
we investigate here using the Community Earth System Model Large Ensemble of climate model simulations. Anthropogenic 
forcing is found to yield large twenty-first-century increases in the frequency of wet extremes, including a more than threefold 
increase in sub-seasonal events comparable to California's 'Great Flood of 1862'. Smaller but statistically robust increases in 
dry extremes are also apparent. As a consequence, a 25% to 100% increase in extreme dry-to-wet precipitation events is pro¬ 
jected, despite only modest changes in mean precipitation. Such hydrological cycle intensification would seriously challenge 
California's existing water storage, conveyance and flood control infrastructure. 


M editerranean climate regimes are renowned for their dis¬ 
tinctively dry summers and relatively wet winters—a glob¬ 
ally unusual combination f Such climates generally occur 
near the poleward fringe of descending air in the subtropics, where 
semi-permanent high-pressure systems bring stable conditions dur¬ 
ing most of the calendar yearl Here, the majority of precipitation 
occurs during the passage of transient storm events during a short 
rainy season^—a distinct seasonality brought about by an equator- 
ward shift in the mid-latitude storm track during winter^. The same 
factors that imbue such regions with their temperate mean climate 
state, however, are also conducive to dramatic swings between 
drought and flood^“^. Subtle year-to-year Jetstream shifts can gen¬ 
erate disproportionately large precipitation variability^—yielding 
highly non-uniform precipitation distributions^ and increasing 
the intrinsic likelihood of hydroclimatic extremes^’^. These effects 
are often amplified in California, where a combination of complex 
topography and over 1,000 km of latitudinal extent yield a great 
diversity of microclimates within the broader ‘dry summer’ regimef 
California’s rapid shift from severe drought to abundant pre¬ 
cipitation (and widespread flooding) during the 2016-2017 winter^^^ 
offers a compelling example of one such transition in a highly popu¬ 
lated, economically critical and biodiverse region^Immediately 
following one of the most intense multi-year droughts on record 
between 2012 and 2016 (refs the state experienced several 
months of heavy precipitation associated with an extraordinarily 
high number of atmospheric river storms during November-March 
2016-2017 (ref. While the heaviest precipitation was concen¬ 
trated in northern Sierra Nevada watersheds, hundreds of roads 
throughout California were damaged by floodwaters and mudslides 
(including a major bridge collapse)In February 2017, heavy run¬ 
off in the Feather River watershed contributed to the failure of the 
Oroville Dam’s primary spillway—culminating in a crisis that forced 
the emergency evacuation of nearly a quarter of a million people^^. 

Previous studies focusing on future changes in California pre¬ 
cipitation have generally reported modest^^"^^ (and/or uncertain) 
changes in regional mean precipitation^’^®. More recent work. 


however, has suggested an increased likelihood of wet years^®"^^ 
and subsequent flood risk^’^^ in California—which is consistent 
with broader theoretical and model-based findings regarding the 
tendency towards increasing precipitation intensity^^ in a warmer 
(and therefore moister) atmosphere^^’^^. Meanwhile, while evi¬ 
dence shows that anthropogenic warming has contributed to 
an increased risk of California drought via increasing tempera- 
tures^^’^^ and increased frequency of seasonally persistent high- 
pressure ridges^’^^’^®’^^ attribution studies focusing directly on 
precipitation have yielded mixed results^^’^^. Contributing addi¬ 
tional uncertainty are climate model simulations suggesting that 
the boundary between mean subtropical drying and mid-latitude 
wetting will probably occur over California^^ potentially yield¬ 
ing strong latitudinal gradients in the precipitation response. 
Thus, while there is already substantial evidence that climate 
change will induce regional hydroclimatic shifts, a cohesive pic¬ 
ture has yet to emerge—presenting serious challenges to deci¬ 
sion-makers responsible for ensuring the resilience of California’s 
water infrastructure^ f 

Importance of large ensemble approach 

We use specific flood and drought events from California’s history 
as baselines for exploring the changing character of precipitation 
extremes. Our use of a large ensemble of climate model simula¬ 
tions^^—the Community Earth System Model Large Ensemble 
(CESM-LENS)—allows us to directly quantify changes in large- 
magnitude extremes. This approach offers a substantial advantage 
over traditional climate model experiments, which yield too small 
a sample size of statistically rare extreme events to draw robust 
inferences without making assumptions regarding the underlying 
precipitation distribution^^ By selecting a wide range of wet, dry 
and dry-to-wet transition (that is, ‘whiplash’) events informed by 
historical analogues, we aim to provide a comprehensive perspec¬ 
tive on the changing risks of regional hydroclimate extremes in a 
manner directly relevant to climate adaptation and infrastructure 
planning efforts. 
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We assess simulated changes in the frequency of California pre¬ 
cipitation extremes caused by increasing atmospheric greenhouse 
gas concentrations. Our overall approach is to (1) determine approx¬ 
imate frequency of occurrence for each event of interest based on 
direct observations or historical accounts, (2) calculate the magni¬ 
tude of events within a preindustrial control climate model simula¬ 
tion that occur with comparable frequency to those observed and 
(3) quantify subsequent frequency changes in similar events under 
a scenario of continued growth in greenhouse gas concentrations. 

We focus on changes in frequency of precipitation events exceed¬ 
ing particular thresholds for two key reasons. First, historical civil 
engineering and risk management practices have been predicated 
on a largely stationary climate^^, and the majority of existing water 
storage and conveyance structures have been constructed under 
such assumptions. Second, a frequency-based approach also offers 
the considerable advantage of implicit climate model bias correc¬ 
tion. While all global climate models exhibit some degree of mean 
precipitation bias in topographically complex California^, the use 
of a long preindustrial control simulation to define return interval 
thresholds allows us to make internally consistent comparisons 
between simulated precipitation distributions at different levels of 
radiative forcing. We can therefore select CESM-LENS precipitation 
thresholds corresponding to approximate return intervals of real- 
world historical events, which serve as analogues for impacts. 

Large increase in extreme wet-event frequency 

We find large, statistically robust increases in the simulated fre¬ 
quency of extremely heavy precipitation events on multiple tim¬ 
escales. All of California experiences a 100-200% increase in the 
occurrence of very high cumulative seasonal precipitation (of a 
magnitude comparable to the 2016-2017 season on a statewide 
basis) by the end of the twenty-first century (Eig. 1). This simu¬ 
lated increase in seasonal wet extremes across California is part 
of a broader regional increase extending across the Pacific coast. 
Seasonal precipitation of this magnitude (equivalent to that asso¬ 
ciated with the 25 year preindustrial control (PIC) return interval) 
has only occurred a handful of times over the lifespan of California’s 
modern water infrastructure. It represents a rarely exceeded but 
not unprecedented threshold, for which there are analogues in the 
recent historical record. 

We note, however, that California’s most severe floods do not 
necessarily coincide with its wettest winters. Instead, regional flood 
events are more directly linked to persistent storm sequences on 
sub-seasonal timescales, which are capable of bringing a signifi¬ 
cant fraction of annual average precipitation over a brief period^’^’^^. 
Thus, to better characterize changes in the frequency of such ‘high 
consequence, low probability’ precipitation events, we use a sub- 
seasonal threshold motivated by the extraordinary sequence of 
‘atmospheric river’ storms that brought extremely severe flooding 
to much of California during the winter of 1861-1862 (refs We 
define 40 day precipitation accumulations exceeding the 200 year 
preindustrial return interval as a measure of occurrences compa¬ 
rable to this benchmark event (see Methods). 

Given the severe impacts even one such occurrence would have 
on California’s existing infrastructure and population centres^^ we 
assess cumulative twenty-first-century risk beginning in the pres¬ 
ent winter season (2017-2018). Figure 2 shows that at least two- 
thirds (66.66%) of LENS ensemble members simulate two or more 
1862-magnitude events over this interval across virtually all of 
California—which represents a dramatic increase in likelihood rela¬ 
tive to preindustrial simulations (where, by definition, the cumu¬ 
lative 83 year likelihood of a single occurrence is less than 50%). 
Strikingly, these findings suggest that California’s major urban cen¬ 
tres (including San Erancisco and Los Angeles) are more likely than 
not to experience at least one such extremely severe storm sequence 
between 2018 and 2060 (Eig. 2b,c) on a business-as-usual emissions 
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Fig. 11 Change in frequency of extremely wet seasons, a, Relative (%) 
change in frequency of extremely wet seasons (meeting or exceeding 
the 25 year PIC return interval for November-March precipitation) at end 
of the twenty-first century (2070-2100, RCP8.5 forcing) relative to the 
preindustrial era (1850 forcing). Cross-hatching signifies 90% statistical 
confidence (P<0.10) in robustness of frequency shifts across the full 
40-member CESM-LENS ensemble. b,c, Time series showing relative (%) 
change in frequency of extremely wet seasons in each year from 1935 to 
2085 (solid green curve) for a cluster of grid boxes in Southern California 
(b) and Northern California (c). Data are smoothed over 30 year intervals, 
and the green shaded range encompasses two-thirds (66.66%) of the 
CESM-LENS ensemble spread (that is, the 16.66th and 83.33th percentile 
bounds). Dashed black horizontal lines in b and c denote zero change in 
frequency. 


trajectory. On a statewide basis, the overall frequency of 1862-mag¬ 
nitude events increases on the order of 300-400% by the end of the 
twenty-first century (Supplementary Eig. 2). 

These increases in extreme wet-event frequency across the large 
ensemble emerge in an essentially monotonic fashion throughout 
most of California beginning between around 2010-2020 for both 
seasonal (25 year) wet events (Supplementary Eig. 3a) and sub-sea¬ 
sonal (200 year) wet events (Supplementary Eig. 3b). In addition to 
the robust ensemble-mean signal by mid-to-late century for both 
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Fig. 2 I Cumulative occurrence of extremely wet sub-seasonal storm 
sequences, a, Cumulative number of extremely wet sub-seasonal storm 
sequences (meeting or exceeding the 200 year PIC return interval for 
cumulative 40 day precipitation) occurring in CESM-LENS between 2018 
and 2100 under RCP8.5 forcing. Cross-hatching signifies regions where at 
least two-thirds (66.66%) of CESM-LENS ensemble members simulate 
two or more such occurrences. b,c, Time series showing cumulative 
number of extremely wet sub-seasonal storm sequences during 2018-2100 
(solid blue curve) for a cluster of grid boxes in Southern California (b) 
and Northern California (c). Blue shaded range encompasses two-thirds 
(66.66%) of the CESM-LENS ensemble spread (that is, the 16.66th 
and 83.33th percentile bounds). Dashed black curve depicts baseline 
cumulative occurrence over an equivalent time interval assuming constant 
preindustrial climate forcings. 


Fig. 3 I Change in frequency of extremely dry seasons, a, Relative (%) 
change in frequency of extremely dry seasons (meeting or falling below 
the 100year PIC return interval for low November-March precipitation) 
at end of the twenty-first century (2070-2100, RCP8.5 forcing) relative 
to the preindustrial era (1850 forcing). Cross-hatching signifies 90% 
statistical confidence (P< 0.10) in robustness of frequency shifts across 
the full 40-member CESM-LENS ensemble. b,c. Time series showing 
relative (%) change in frequency of extremely dry seasons in each year 
1935-2085 (solid brown curve) for a cluster of grid boxes in Southern 
California (b) and Northern California (c). Data are smoothed over 
30 year intervals, and the brown shaded range encompasses two-thirds 
(66.66%) of the CESM-LENS ensemble spread (that is, the 16.66th and 
83.33th percentile bounds). Dashed black horizontal lines in b and c 
denote zero change in frequency. 


extreme wet events, these results also imply that an increased likeli¬ 
hood of large precipitation accumulations relative to a preindustrial 
climate may already exist. However, we note that a handful of out¬ 
lying ensemble members suggest a chance that internal variability 
could delay emergence of an attributable signal by several decades 
(Fig. lb,c and Fig. 2b,c). 

Shifts in extreme dry-event frequency 

The simulated frequency of extremely dry years also increases across 
nearly all of California (Fig. 3). An increased likelihood of these 


extremely dry rainy seasons (that is, exceeding the 100 year return 
interval, analogous to the 1976-1977 drought and slightly drier than 
2013-2014) first emerges weakly across portions of the state as early 
as the 1980s, before emerging in a statistically robust manner across 
Northern California around 2010-2020 and Southern California 
later in the century (around 2060; Supplementary Fig. 4a). Notable 
are especially large increases (>140%) in frequency that occur 
across Southern California after 2050 (Fig. 3b and Supplementary 
Fig. 4a), though we emphasize that substantial increases on the 
order of -1-80% extend across most of Northern California. Except 
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for southernmost California, much of this increased risk emerges 
during the twentieth and early twenty-first centuries (Fig. 3b,c and 
Supplementary Fig. 4a)—suggesting that the likelihood of individ¬ 
ual dry seasons may already be increased relative to the preindus¬ 
trial period. 

In contrast, changes in the occurrence of extremely dry consecu¬ 
tive years analogous to the record-low 3 year cumulative statewide 
precipitation observed during 2013-2015 (100 year return interval 
on a 3 year basis) exhibit a more complex temporal and spatial struc¬ 
ture (Supplementary Fig. 4b). By the end of the twenty-first cen¬ 
tury, only far southern California experiences a robust increase in 
the frequency of consecutively dry seasons (Supplementary Fig. 2), 
while the rest of California does not experience statistically signifi¬ 
cant changes. Further analysis shows that this divergence between 
single and consecutive dry-season frequency shifts arises from 
the increased pace of future wet-year increases relative to dry- 
year increases, which is especially pronounced across Northern 
California (Supplementary Fig. 5c). These findings suggest that 
future multi-year droughts in California may exhibit an increased 
propensity to be interrupted by very wet interludes. 

Emergence of 'precipitation whiplash' signal 

Given the large simulated increase in the frequency of both dry 
and wet extremes, we test whether the frequency of rapid transi¬ 
tions between dry and wet conditions—similar to the precipitation 
whiplash that occurred between the recent 2012-2016 drought and 
2016-2017 floods—also increases. For the purposes of this study, we 
define precipitation whiplash as the occurrence of two consecutive 
years during which rainy season (November-March) precipitation 
falls under the PIC 20th percentile (in the first year) and subse¬ 
quently exceeds the PIC 80th percentile (in the following year). 
Figure 4 confirms such an increase during the twenty-first century 
throughout California. We report a strong latitudinal gradient in 
the year-to-year (interannual) whiplash response to anthropogenic 
forcing, ranging from an -100% increase across Southern California 
to an -25% increase in Northern California (Fig. 4a). These whip¬ 
lash increases first emerge in the south relatively early in the twenty- 
first century (2010-2020) before spreading progressively northward 
in a statistically robust manner in the following decades (especially 
after 2050; Supplementary Fig. 5). We also investigate changes in 
month-to-month (sub-seasonal) precipitation variability during the 
canonical wet season. We report modest but widespread increases 
of 20-30% across a broad swath of the northeastern Pacific region, 
again extending across all of California (Supplementary Fig. 2c). 

Together, these shifts represent a marked increase in both the 
interannual and intraseasonal variability of precipitation, especially 
in Southern California—which is noteworthy for two distinct rea¬ 
sons. First, natural precipitation variability in this region is already 
large^, and projected future whiplash increases would amplify exist¬ 
ing swings between dry and wet years (and between dry and wet 
months within the rainy season). Second, the robust emergence of 
a precipitation whiplash signal across a wide range of timescales 
(Fig. 4 and Supplementary Fig. 2c) is remarkable, as large-scale 
atmospheric variability over the North Pacific relevant to California 
precipitation is dominated by different physical processes and 
associated remote teleconnections on sub-seasonal (for example, 
the Madden-Julian Oscillation^^) and interannual (for example, El 
Nino-Southern Oscillation (ENSO)^^) timescales. 

Increasingly sharp seasonality of California wet season 

We find a distinct sharpening of Californias future mean sea¬ 
sonal cycle (Fig. 5). While winter mean precipitation increases 
across most of California, mean precipitation during autumn 
(September-November) and especially spring (March-May) 
decreases nearly everywhere, which is consistent with previous 
findings from the Coupled Model Intercomparison Project Phase 
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Fig. 4 I Change in frequency of precipitation whiplash events, a, Relative 
(%) change in frequency of dry-to-wet precipitation whiplash events 
(years with November-March precipitation at or above the 80th 
preindustrial percentile immediately followed by a year with precipitation 
below the 20th preindustrial percentile) at end of the twenty-first century 
(2070-2100, RCP8.5 forcing) relative to the preindustrial era (1850 forcing). 
Cross-hatching signifies 90% statistical confidence (P< 0.10) in robustness 
of frequency shifts across the full 40-member CESM-LENS ensemble. 
b,c, Time series showing relative (%) change in frequency of whiplash 
events in each year 1935-2085 (solid purple curve) for a cluster of grid 
boxes in Southern California (b) and Northern California (c). Data are 
smoothed over 30 year intervals, and the purple shaded range 
encompasses two-thirds (66.66%) of the CESM-LENS ensemble spread 
(that is, the 16.66th and 83.33th percentile bounds). Dashed black 
horizontal lines in b and c denote zero change in frequency. 


5 (CMIP5) ensemble^^ This striking contrast between the drying 
marginal and wetting core rainy season months results in a large 
ensemble-mean increase (35% to 85% from north to south) in the 
ratio of overall wet season precipitation falling between November 
and March relative to cumulative precipitation during the four 
months of September, October, April and May (Fig. 5). This 
increase in sharpness of precipitation seasonality suggests that the 
already distinct contrast between Californias long, dry summers 
and relatively brief, wet winters will probably become even more 
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Fig. 5 I Shifts in precipitation seasonality. Relative changes in CESM- 
LENS monthly mean precipitation at the end of the twenty-first century 
(2070-2100) as a percent of the PIC climatology for each calendar 
month for a range of latitudes spanning the California coast. Percentages 
in the legend denote relative changes in mean 'seasonal sharpness' at 
each latitude, defined as the ratio between precipitation falling during 
the core rainy season (November-March; blue background shading) to 
that cumulatively falling during the marginal rainy season (September- 
October, April-May; red background shading). Curves are colour coded by 
latitude (and therefore by mean seasonal precipitation, which increases 
monotonically with latitude). Dashed black horizontal line denotes zero 
change in magnitude. 

pronounced during the twenty-first century While a comprehen¬ 
sive assessment is beyond the scope of this study, we note that 
autumn and spring drying trends have recently begun to emerge in 
observations across California (Supplementary Fig. 7)—suggest¬ 
ing that the projected concentration of precipitation into an even 
narrower season may already be underway 

Increase in extremes despite modest mean change 

The substantial increases in California precipitation extremes over 
a wide range of timescales and intensities occur despite only modest 
changes in rainy season mean precipitation. By the end of the twenty- 
first century, the CESM-LENS ensemble mean depicts modest cool- 
season wetting over the northern portion of the state (<20-30%; 
Fig. 5), with little change in the south (~0%; Supplementary Fig. 6)— 
similar to the CMIP5 multi-model ensemble-mean response 
(Supplementary Fig. 6). Yet over the same interval, the frequencies 
of both extreme dry seasons and whiplash events increase by over 
50% over much of the state (Figs. 3 and 4), and the frequency of 
extreme wet events increases by well over 100% nearly everywhere 
(Fig. 1 and Supplementary Fig. 2). This remarkable divergence 
between simulated future mean and extreme climate is especially 
pronounced across Southern California. For example: simulated 
mean precipitation in Southern California exhibits little trend by the 
end of the twenty-first century (Supplementary Fig. 5), despite an 
-200% increase in extremely dry seasons (Fig. 3^ an -150% increase 
in extremely wet seasons (Fig. 1), a >500% increase in extreme sub- 
seasonal precipitation events (Fig. 2 and Supplementary Fig. 2) and 
an -75% increase in year-to-year whiplash (Fig. 4). Importantly, 
these findings suggest that lingering uncertainty regarding the mag¬ 
nitude (and even sign) of regional mean precipitation change does 
not preclude statistically robust conclusions regarding changes in 
precipitation extremes. 

Changes in processes responsible for extremes 

We confirm that simulated large-scale atmospheric circulation 
patterns associated with California wet and dry extremes are 


substantially similar to those observed historically (Supplementary 
Fig. 12). Wet years are linked to strong low-pressure anomalies 
over the northeastern Pacific Ocean (Fig. 6a), with downstream 
enhancement of the storm track just west of California^’^^. Dry 
years coincide with seasonally persistent high pressure extend¬ 
ing across the northeastern Pacific (Fig. 6d), which reinforces the 
climatological mean ridge along the West Coast^^ and prevents 
storms from reaching California^^. 

We find that future (representative concentration pathway 8.5 
(RCP8.5)) wet and dry years are linked to broadly similar atmospheric 
circulation anomalies as those in the past (PIC; Fig. 6b,e)— suggest¬ 
ing that the spatial character of large-scale features driving California 
precipitation extremes may remain relatively stationary. However, 
given that subtle storm track perturbations can yield disproportion¬ 
ately large shifts in California precipitation^’^^’^^^, we note two poten¬ 
tially important differences between the RCP8.5 and PIC composites. 
In future wet years, low pressure over the North Pacific is deeper to 
the west of California (Fig. 6b), which previous work using models 
in the CMIP5 experiment has suggested is linked to a localized east¬ 
ward extension of the jetstream^’^^^. In RCP8.5 dry years, atmospheric 
pressure anomalies in the immediate vicinity of California are similar 
to PIC dry years (Fig. 6e), but are higher in adjacent regions—sug¬ 
gesting a broader, more longitudinally oriented atmospheric ridge 
pattern and subsequent poleward storm track shift. 

We also report large increases in atmospheric water vapour dur¬ 
ing both future wet and dry years. While this moistening is not in 
itself surprising—given the well-understood thermodynamic con¬ 
sequences of the Clausius-Clapeyron relation^^—we point out that 
RCP8.5 dry years occur in an atmosphere moister than that during 
even the wettest years of the PIC simulation (Fig. 6c,f). The fact that 
California dry years occur more frequently (Fig. 3) suggests that 
simulated (thermodynamic) moistening must be counteracted (at 
least periodically) by changes in the frequency and/or intensity of 
atmospheric circulation patterns that prevent precipitation-bearing 
storms from reaching California, such as persistent high-pressure 
systems^^ or transient poleward shifts in the East Pacific storm tracks 
Conversely, the (thermodynamic) increase in water vapour would 
probably reinforce the (dynamic) effect of deeper North Pacific low 
pressure during wet years—which may underlie the relatively larger 
simulated increase in extreme wet-event frequency (versus dry- 
event frequency). Nonetheless, we emphasize that further work is 
necessary to better understand underlying changes in both remote 
(tropicaP^’^^ and Arctic teleconnections^^“^^) and regional-scale (that 
is, atmospheric rivers^ and orographic precipitation) influences. 

Societal implications of hydrological intensification 

Collectively, our findings suggest that anthropogenic warming will 
bring about large increases in the frequency of California hydrocli- 
matic extremes similar or greater in magnitude to those that have 
historically caused widespread disruption. These changes in the 
character of California precipitation emerge in a large single-model 
ensemble despite only modest trends in mean precipitation— 
strongly suggesting that the regions already variable year-to-year 
climate is likely to become even more volatile. Historically observed 
impacts of droughts and floods may in many cases offer reason¬ 
able analogues for the human and environmental impacts of future 
events of a similar magnitude, but increasingly wide swings between 
dry and wet conditions will threaten to upset the already precarious 
balance between competing flood control and water storage imper¬ 
atives in California. 

Moreover, we report a substantial increase in the projected risk 
of extreme precipitation events exceeding any that have occurred 
over the past century—meaning that such events would be unprece¬ 
dented in California’s modern era of extensive water infrastructure. 
Eew of the dams, levees and canals that currently protect millions 
living in California’s flood plains and facilitate the movement of 
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Fig. 6 I Large-scale atmospheric conditions linked to California precipitation extremes. a,d, Composite maps depicting anomalies in November-March 
500 hPa GPH during wet (80th percentile, left column) and dry (20th percentile, right column) March-November seasons in the CESM-LENS simulation. 
b,e, RCP8.5 (2081-2100) composite anomaly patterns for GPH (see Methods). c,f. Difference between RCP8.5 (2081-2100) and PIC composite anomaly 
patterns for column-integrated water vapour. 


water from Sierra Nevada watersheds to coastal cities have been 
tested by a deluge as severe as the extraordinary 1861-1862 storm 
sequence—a repeat of which would probably lead to considerable 
loss of life and economic damages approaching a trillion dollars^^’^®. 
Our results suggest that such an event is more likely than not to 
occur at least once between 2018 and 2060, and that multiple occur¬ 
rences are plausible by 2100 on a business-as-usual emissions trajec¬ 
tory Therefore, recognizing that risks associated with hydroclimatic 
extremes may rise more rapidly than the gradual projected shift in 
regional mean precipitation might otherwise suggest will be a criti¬ 
cal step in ensuring resilience amid a warming climate. 

Methods 

Methods, including statements of data availability and any asso¬ 
ciated accession codes and references, are available at https://doi. 
org/10.1038/s41558-018-0140-y 
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Methods 

Datasets used in this study. The CESM-LENS is a large ensemble of fully coupled 
model simulations designed to explore multiple realizations of internal climate 
system variability on long timescales. We used precipitation output from an 
1,800 year PIC run and 40 separate simulations of the twentieth century (20C; 
1920-2005) and RCP8.5 (2005-2100) climate change scenario^^ Each of the 40 
20C + RCP8.5 realizations is generated using the same climate model but with 
slightly perturbed initial conditions, which yield different time sequences of daily- 
to decadal-scale internal variability"^”. Thus, CESM-LENS offers an opportunity 
to examine robust changes in extreme events across a wide range of simulated 
internal variability—a considerable advantage relative to other investigations that 
have historically been limited by the shortness of the observational record and the 
relative infrequency of extreme hydroclimatic events in smaller ensembles. 

We also used precipitation output from climate model simulations generated 
as part of the CMIP5 projecff*^ for comparison with CESM-LENS simulations. We 
constructed a multi-model ensemble consisting of 78 realizations from 35 distinct 
climate models, where each distinct model receives equal weight in the ensemble 
mean and fields are interpolated to a common 2.5° grid. 

We used the National Climatic Data Center’s nClimDiv observational 
divisional dataset to determine the relative rank of historical precipitation events 
to estimate approximate return intervals (Supplementary Eig. 2). Existing biases 
between simulated and observed precipitation were implicitly accounted for 
using the methodology described below. Einally, we used gridded observational 
precipitation data to validate CESM-LENS precipitation (see section on Suitability 
of CESM-LENS for simulating California precipitation). 

Quantifying changes in frequency of extreme hydroclimatic events. We seek 
to quantify changes in the frequency of extreme wet, dry and whiplash events 
in a manner that implicitly accounts for model biases and sidesteps parametric 
assumptions regarding the underlying shape of the precipitation distribution. In 
doing so, we focus on relative changes in the frequency of occurrence of events 
exceeding fixed thresholds defined using the PIC simulation. PIC atmospheric 
greenhouse gas concentrations are maintained at constant levels similar to those 
before the start of the Industrial Revolution (that is, year 1850 levels)—representing 
a counterfactual climate without human influence. We subsequently compare the 
relative change in frequency of specific events in the 20C-I-RCP8.5 simulations. 

The 20C forcing includes rising greenhouse gas and aerosol concentrations close 
to those observed in the historical record, and RCP8.5 forcing includes projected 
greenhouse gas increases between 2005 and 2100 based on a business-as-usual 
emissions trajectory^"^. We examine the RCP8.5 simulations as they most closely 
resemble the observed emissions trajectory to date"^T and they provide a larger 
signal-to-noise ratio for statistically rare extreme events^^. 

We restrict our analysis of extremes to the months of November to March, 
representing the peak of the California rainy season^^. While the seasonal peak 
of monthly precipitation occurs earlier in Northern California than in Southern 
California, extreme winter-like precipitation events associated with mid-latitude 
cyclonic activity can occur during any of these calendar months throughout the 
state. Thus, data for November-March are pooled to create a single, spatially 
explicit rainy season distribution of precipitation accumulations in each grid box. 

Using direct observational (and indirect historical) records, we first estimate 
the approximate return interval for each event of interest. Our definition of return 
interval (R) is consistent with that widely used in the climate and civil engineering 
literature^”: the likelihood of occurrence in any given year of an event with an 
nyear return interval is H n. Using this fixed return interval, we then calculate 
the precipitation value (p) in the CESM-LENS PIC simulation that occurs with 
frequency 1 / n at each climate model grid box. In the 40-member ensemble, we 
then count exceedances of the specific return values and normalize the PIC and 
20C -I- RCP8.5 runs on an events-per-year basis. 

Once these counts are tabulated across all ensemble members and each season, 
we calculate percent changes in the frequency F of a given event at time t relative to 
the PIC period: 

4F,= & (1) 

^PIC 

This approach is modelled after ref. and references therein. We note that our 
use of a fixed return interval, rather than percentile values, allows us to directly 
compare precipitation extremes spanning timescales from monthly to multi¬ 
annual. Eor example: the 98th percentile of daily precipitation might be expected to 
occur several days per year, but the 98th percentile of annual precipitation might be 
expected to occur only twice per century. 

Definition of precipitation whiplash. Given the potential for adverse human 
impacts of rapid transitions between dry and wet conditions in California (as 
occurred during 2016-2017), we formalize two precipitation whiplash metrics 
across a range of timescales. We define year-to-year (interannual) whiplash years 
to be those during which seasonal (November-March) precipitation accumulation 
meets or exceeds the PIC 80th percentile and which were immediately preceded by 
a year with seasonal precipitation at or below the PIC 20th percentile. We define 


within-season (intraseasonal) whiplash as the standard deviation of monthly 
precipitation within individual rainy seasons across the ensemble. 

Selection of extreme events and climate model analogues. 2016-2017 (wet) 
analogue with 25year return interval. This threshold is based on the extremely 
wet 2016-2017 winter, during which record wet seasonal precipitation (return 
interval >100 years) occurred across portions of Northern California and relatively 
smaller positive anomalies occurred over Southern California (return interval 
-lOyears), yielding a mean statewide return interval of approximately 15-25yr for 
November-March. Given the strong latitudinal gradient in relative abnormality of 
seasonal precipitation during this period—and the consequent variation in adverse 
societal impacts across the state—we use a 25 yr return interval as a compromise 
threshold. 

1861-1862 (wet) analogue with 200year return interval. This threshold is based 
on the extraordinary sequence of atmospheric river storm events that brought 
extremely severe and widespread flooding to much of California during the 
1861-1862 winter. Much of what is known about the ‘Great Elood of 1862’ has 
been pieced together from informal historical accounts and newspaper records 
from the then-nascent State of California^®’”^. Such records suggest that the most 
intense period of nearly continuous precipitation occurred between late December 
1861 and late January 1862 over an approximately 40 day period, yielding rainfall 
accumulations over 1 m in some locations^^ While this event occurred before 
the advent of reliable meteorological observations in California, palaeoclimate 
evidence from sediment records in coastal river systems suggests that events 
comparable to the 1861-1862 flood are associated with an approximately 200 year 
return interval”. Given the substantial uncertainties regarding the exact duration, 
magnitude and recurrence interval of the 1861-1862 event, we define our analogue 
as the 40 day cumulative precipitation during all November-March periods with a 
return interval of 200 years in the PIC simulation. The magnitude of this event is 
larger than that envisioned in the ‘ARkStorm’ natural hazard contingency planning 
scenario jointly developed by the United States Geological Survey and the State of 
California^® but smaller in magnitude compared with several other probable events 
in the past millennium”. 

Recent research has suggested that a modern recurrence of the 1861-1862 
flood would probably have a catastrophic human and socioeconomic tolP®’^”. 

Thus, despite the fact that such an event remains unlikely in any given year even 
under strong greenhouse forcing, a multi-fold relative increase in physical event 
likelihood combined with a high degree of socioeconomic vulnerability collectively 
yield a substantial increase in the overall risk associated with such an event over a 
period of decades”. 

1976-1977 (dry) analogue with 100year return interval. This threshold is based on 
the extremely dry conditions that occurred during winter 1976-1977, which was 
the driest such period in California’s 122 year observational record. This short-lived 
but intense drought led to acute water shortages in regions dependent on surface 
runoff from smaller hydrological basins and without direct access to State Water 
Project or Central Valley Project reservoirs. We conservatively assume a 100year 
return period for this single-year event analogue threshold. 

2012-2016 (dry) analogue with 100year return interval. This threshold is based 
on the multi-year drought that occurred in California between late 2012 and early 
2016. A substantial fraction of overall drought magnitude and associated impacts 
can be attributed to extremely warm temperatures that coincided with the lack of 
precipitation during successive winters^*^ -®, and while 2013 was the driest calendar 
year on record in California®, no individual November-March period was the 
drier than 1976-1977. Nevertheless, the driest consecutive 3 year period (and 
consecutive November-March seasons) on record in California did indeed occur 
between 2013 and 2016 (ref. ”), and we use this 3-year threshold as a benchmark 
for a high-impact, multi-year drought. We emphasize that the widespread 
environmental impacts of this event were substantially exacerbated by record 
warm temperatures—which are expected to be a signature of future droughts in 
this region as the climate warms"® ””. Consistent with the single-year dry event, we 
conservatively assume a 100 year return interval. 

Quantification of statistical significance. All figures showing spatial changes 
(latitude-longitude and time-latitude maps) represent 30 year running means in 
the 20C-I-RCP8.5 simulations. In all significance assessments below, these 30year 
mean changes in event frequency are compared with resampled (bootstrapped) 
time series from the PIC simulation. Confidence intervals represent climate change 
signals that fall outside the sampled range of PIC variability with 90% confidence, 
representing a high statistical bar given the very wide range of simulated internal 
variability that exists within CESM-LENS””. 

Change-in-frequency maps. To provide a robust measure of statistical significance 
for change-in-frequency maps (Eigs. 1-4 and Supplementary Eig. 2), we used a 
bootstrap resampling approach. Eor seasonal extreme events (25 year wet events 
and 100year dry events), we generated 10,000 random time series by selecting wet 
season precipitation totals from the full 1,800year PIC run (with replacement). 
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The length of each time series corresponds to the return interval of the event 
itself: in this instance, either 25 or 100 years. Next, in each resampled time series, 
the number of exceedances of 25 or 100 year events was calculated and translated 
into a ratio relative to the full PIC count, producing a distribution of 10,000 
bootstrapped ratio values at each grid point. Finally, these distributions were used 
to determine the rarity of the simulated ratios relative to internal variability in the 
PIC run. In all map plots, cross-hatching for wet (dry) events represent a value at 
or above (at or below) the 90th (10th) PIC percentile—signifying 90% confidence 
(P< 0.1) that ratios fall outside of the PIC internal variability. In time series plots 
shown in Figs. 1-4, error bars represent two-thirds of the spread (±33.3%) of 40 
ensemble members, and the average is calculated on a 30 year running mean basis 
to distinguish long-term trends from interannual variability. 

We used a modified version of the bootstrap resampling approach described 
above for the 200 year event change-in-frequency map (Fig. 2 and Supplementary 
Fig. 2). Given the computational constraint of large-N resampling using this 
dataset, 100 resampled time series were constructed for this particularly extreme 
event (we confirm that N= 100 is a sufficiently large sample size for estimates 
of precipitation values to stabilize). For each time series, 200 years were chosen 
at random (each of which contains a distribution of 40 day running sums 
of November-March precipitation). We note that temporal autocorrelation 
can become problematic when counting occurrences using 40 day running 
sums for high-magnitude events. To ensure that our frequency counts do not 
unintentionally count these extremes twice, our algorithm skips ahead by 40 days 
each time a 200-year-magnitude event is encountered. As for other extreme events, 
we calculated the ratio of event frequencies between the full PIC run and the 
bootstrapped time series at each grid point. 

For change-in-frequency maps of seasonal whiplash events (Fig. 4a), 10,000 
block bootstrapped time series were generated using randomly chosen segments of 
100 consecutive years. This consecutive-year (block) approach is necessary because 
our whiplash definition depends on sequential dry-to-wet transitions; thus, to 
appropriately sample the internal variability, the underlying temporal sequence in 
the PIC simulation must be preserved. 

For change-in-frequency maps of month-to-month precipitation variability 
(Supplementary Fig. 2c), P values are calculated using 10,000 bootstrapped 
40 year time series of November-March PIC precipitation. In each bootstrap 
iteration, 40 years were selected at random and with replacement, and the standard 
deviation across 200 (40 X 5) November-March model months was calculated 
for comparison with the CESM-LENS 40-member ensemble. Distributions of 
these month-to-month variability measures were generated at each grid point and 
normalized by the full PIC run to represent ratios. 

For time series plots referring to either Southern California or Northern 
California, we use a spatially smoothed mean value (defined as the average value 
within a 3 X 3 grid box cluster centred on the original CESM-LENS grid box closest 
to the actual latitude/longitude of Los Angeles and San Francisco, respectively). 

Time-latitude plots. The bootstrapped time series discussed above are also used 
to calculate statistical significance for time-latitude plots. After these bootstrapped 
time series were generated, further calculations were completed using the average 
of three contiguous west-to-east-oriented grid boxes along the California coast 
at each relevant latitude. As before, P values <0.1 resulting from the significance 
test imply rejection of the null hypothesis that the ratios in the time-latitude plots 
are within the range of simulated internal variability in the PIC simulation. In 
the present case, rejection of the null hypothesis is interpreted to mean that the 
20C + RCP8.5 distribution is statistically distinguishable from the climate without 
humans’ control. 

Analysis of large-scale atmospheric conditions linked to extremes. We created 
anomaly composite maps for 500 hPa geopotential heights (GPH; Fig. 6a,d) and 
column-integrated water vapour (Fig. 6c,f) for wet years (exceeding the 80th 
preindustrial percentile) and dry years (falling below the 20th percentile) in the 
LENS simulations during both the full 1,800 year PIC simulation and across the 
40 ensemble members of the RCP8.5 simulation between 2081 and 2100 (yielding 
a sample size of 760 model years). To simplify visual comparison of the anomaly 
patterns (where spatial gradients determine the geostrophic wind field) a spatially 
constant component of the thermal dilation^’^'^ of the atmosphere (Supplementary 
Fig. 8a) was removed by subtracting the field mean difference in GPH (RCP8.5 
minus PIC for a broad region near California (20° N-60° N, 180° E-250° E) 
from all grid points before generating RCP8.5 anomaly fields (Eig. 6a,d). We also 
compare the difference in anomalies between the RCP8.5 and PIC for wet and 
dry years (Pig. 6b,c,e,f) to emphasize the relative similarity of the underlying 
atmospheric circulation features between these periods. 

Suitability of CESM-LENS for simulating California precipitation. The majority 
of findings in this study were derived using a global climate model operating 
at relatively coarse spatial resolutions (nominally 1°, or -100 km, for the CESM 
configuration in LENS^^). Given the importance of fine-scale topography in 
influencing California precipitation extremes^ the relative spatial coarseness of 
model data used in this study precludes quantitative estimates of future runoff 
volume and flood flows at the watershed scale. Indeed, we note that that a local 
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minimum in the relative increase in wet extremes exists along the axis of maximum 
topographic slope in the CESM-LENS boundary conditions (Supplementary 
Pig. 9b), suggesting a possible nonlinearity in orographic precipitation scaling with 
warming (a possibility supported by recent high-resolution modelling experiments, 
for example, ref. ^^). In the present study, however, our focus on relative (rather 
than absolute) changes in the frequency of various precipitation extremes across 
broad regions implicitly accounts for possible simulated precipitation biases arising 
from coarse model resolution and other sources. 

We reiterate that CESM-LENS, while state-of-the-art, is a single-model 
ensemble. Substantial intermodel differences do exist in the simulated atmospheric 
response to anthropogenic forcing/ but we have chosen to focus on results using 
exclusively CESM-LENS data for two reasons. Pirst, CESM-LENS includes 40 
ensemble members over a 180year simulation (1920-2100)—yielding a very 
large (7,200-model year) sample size, allowing us to directly examine very rare 
extreme events (such as the ‘200 year flood’) without making assumptions about 
the character of the underlying statistical distribution. Combining this with the 
1,800 year PIC simulation allows for evaluation of very rare events that would 
otherwise not be possible in the observational record or CMIP5 experiments. 
Recent evidence suggests that (1) CESM reproduces both the mean and variability 
of observed California precipitation with reasonable fidelity‘s (Supplementary 
Pig. 10), (2) the CESM-LENS single-model ensemble mean lies close to the 
median of the CMIP5 multi-model ensemble mean in the vicinity of California 
(Supplementary Pig. 6) and (3) CESM reproduces remote teleconnections 
(that is, those associated with ENSO) critical to California precipitation^^ 
(Supplementary Pig. 11). 

To independently confirm that CESM-LENS is an appropriate tool for 
investigating changes in California precipitation, we compare simulated versus 
observed precipitation distributions for the 20C historical period for Northern and 
Southern California. We perform this validation using a gridded observational 
dataset (the Global Precipitation Climatology Project (GPCP) version 2.3^®). To 
generate distributions for Northern and Southern California regions (depicted 
as white boxes in Supplementary Pig. 9), we sum November-March precipitation 
over land-only grid boxes with centroids that fall within a 200 km radius of San 
Prancisco and Los Angeles, respectively, during the 1980-2016 period of mutual 
overlap between the CESM-LENS 20C simulations and the GPCP dataset. 

We find that despite slight positive bias in median seasonal precipitation (<5% 
for Southern California and+9% for Northern California), the overall shape of the 
CESM-LENS distribution for both regions is statistically indistinguishable from 
observations at the 5% level using a Kolgoromov-Smirnoff test (Supplementary 
Pig. 10). We note that the CESM-LENS distribution tends to have slightly longer 
tails than observations, but this is unsurprising given that the effective sample 
size of the historical simulation (1,440 model years) is much larger than for 
the observational dataset (40 years) and the range of distribution among the 40 
realizations encompasses the GPCP data (see horizontal bars in Supplementary 
Pig. 10). As our overall methodology implicitly accounts for mean biases in 
precipitation, and our focus is on extreme events in the upper and lower the 
tails of the distribution, the outcome of this validation exercise strongly suggests 
CESM-LENS is capable of capturing both the median and underlying interannual 
variability of California hydroclimate. 

Having confirmed these measures of fidelity of simulated California 
precipitation in CESM-LENS relative to observations, we assess whether the 
ensemble also captures the large-scale physical processes and teleconnections 
responsible for precipitation variability in this region. To test whether CESM- 
LENS plausibly reproduces the observed ENSO teleconnection, we perform linear 
regression of 500 hPa GPH on sea surface temperatures in the ‘Nino3.4’ region of 
the tropical eastern Pacific Ocean using data from both CESM-LENS and National 
Centers for Environmental Prediction/National Center for Atmospheric Research 
(NCEP/NCAR) RH*^. The spatial pattern of the GPH teleconnection is substantially 
similar between CESM-LENS simulations and the R1 reanalysis (corroborating 
results previously shown in ref. and is characterized by a deepening of North 
Pacific low pressure and a more modest decrease in mid-tropospheric heights 
eastward over California and the southern tier of the United States (Supplementary 
Pig. 1 la,b) during El Nino events. We further confirm that the mean position and 
magnitude of the cool-season (November-March) jetstream is close to that in 
observations (Supplementary Pig. 1 lb,c), though we also point out that even subtle 
biases could potentially lead to latitudinal shifts in the location of precipitation 
extreme changes discussed here (an issue that has been raised in previous work^°). 
Pinally, the large-scale atmospheric circulation patterns during California wet and 
dry years, respectively, exhibit similar spatial patterns and magnitudes to those 
observed during historical wet and dry years (Supplementary Pig. 12), especially 
in key regions near the US West Coast. Collectively, these results suggest that 
CESM-LENS is an appropriate tool for use in characterizing changes in regional 
precipitation extremes in the vicinity of California. 

Code availability. The code used in the analyses described in this study is available 
from the corresponding author upon reasonable request. 

Data availability. Precipitation data from the CESM-LENS simulations are 
available from the University Corporation for Atmospheric Research (http://www. 
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cesm.ucar.edu/projects/community-projects/LENS/data-sets.html). Precipitation 
data for California are available from the National Oceanic and Atmospheric 
Administration National Climatic Data Center (nClimDiv, www.ncdc.noaa.gov/ 
monitoring-references/maps/us-climate-divisions.php) and National Oceanic and 
Atmospheric Administration Earth System Research Laboratory (NOAA ESRL) 
(GPCP, https://www.esrl.noaa.gov/psd/data/gridded/data.gpcp.html). CMIP5 
ensemble data were obtained from Lawrence Livermore National Laboratory’s 
Earth System Grid portal (https://esgf.llnl.gov) via the Royal Netherlands 
Meteorological Institute Climate Explorer (https://climexp.knmi.nl). Geopotential 
height and wind data from NCEP/NCAR R1 are available from NOAA ESRL 
(https://www.esrl.noaa.gov/psd/data/gridded/data.ncep.reanalysis.html) and 
additional composite data were created using the NOAA ESRL plotting tool 
(https://www.esrl.noaa.gov/psd/ cgi-bin/data/composites/printpage.pl) . 
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